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WATSON RESUMMATION OF A CLASS OF 
HAUSDORFF TRANSFORMED POWER SERIES 

E. DE MICHELI AND G. A. VIANO 



Abstract. In this paper we study a class of Hausdorff— transformed power 
series whose convergence is extremely slow for large values of the argument. 
We perform a Watson-type resummation of these expansions, and obtain, by 
the use of the Pollaczek polynomials, a new representation whose convergence 
is much faster. We can thus propose a new algorithm for the numerical eval- 
uation of these expansions, which include series playing a relevant role in the 
computation of the partition function in statistical mechanics. By the same 
procedure we obtain also a solution of the classical Hausdorff moment problem. 



1. Introduction 



Following Fuchs, Rogosinski [HJ^] and Hardy [5] we say that f(x) is summable 
to the value S by the continuous Hausdorff method (or f(x) is IK-summable to S) 

i2 : ff(*o= / f(xt)d X (t) (x>o) (i) 

m ■ J ° 

tends to S as x — > oo. Integral Qj is a Lebesgue-Stieltjes integral, f(x) is Borel 
summable and bounded in every finite interval [0, x], and x(t) is of bounded varia- 
tion in [0, 1]. Next, in Hardy 2J the following theorem is proved. 

O ■ 

^5 Theorem 1 (Hardy). In order that the transformation Q should be regular, i.e., 

that f(x) — ► S (for x — » oo) should imply g(x) — > S (for x — * oo), it is necessary 
and sufficient that x(0 + ) = x(0) = and x(l) = 1. 

H ! 

Let us now suppose that f(x) is a function of x regular on the positive real axis, 
^ ■ and so expressible in the form 

f(x) = ±^x^±a n x-. (2) 

n=0 ' ra=0 

Then, we substitute expansion J2} in integral Q and, in view of the uniform con- 
vergence of the Taylor series, we can exchange the sum with the integral and obtain 



J oo oo „i 

g(x) = / T a n (xt) n d X (t) = V a n x n 

JO n „-n JO 



t n d X (t). (3) 



The last integral at the r.h.s. of © represents the Hausdorff moment fj, n , i.e., 

Mn= / t n d X (t)= f t n u(t)dt. (4) 
Jo Jo 

In formula x(i) is supposed to be a real function of bounded variation in 
t € [0, 1], and the numbers pL n are called moment constant, of rank n, of \- If wc 

l 
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suppose, without loss of generality, that x(0) = 0, x(l) = 1, and x(0 + ) = x(0) = 0, 
so that x(~t) is continuous at the origin, then /x„ is called a regular moment constant 
(see Theorem Q] above). Moreover, the following theorem can be proved UJ. 

Theorem 2 (Hardy). Sums, differences and products of moment constants are 
themselves moment constants. The product of two regular moment constants is a 
regular moment constant. 

Two relevant examples of regular Hausdorff transformations are the following 9, 
Theorem 200]: 

(i) p tn =lj\ n {l-t)^- 1 Ut= + ' (£>0), (5) 

which corresponds to the Cesaro transformation ; 

which corresponds to the Holder transformation flw. 

From J2J) we are naturally led to consider expansions of the following form: 

£^W", P) 

where the terms \i n are Hausdorff moments. If we suppose that / (n) (0) = (-1)", 
we obtain expansions which read: 

9{x) = Y,-^-^ n - (8) 

n=0 

These series may be slowly convergent, for values of x sufficiently large. We thus 
face a serious problem of numerical analysis, which is quite relevant in view of 
the fact that sums like expansion JSJl occur in several problems, including some of 
physical interest. For instance, 

(a) The Laplace transform of the functions of compact support gives rise to 
sums of the form (JSJ, if we expand in series the exponential e~ xt . This 
case appears in statistical mechanics, where the partition function is the 
Laplace transform of the density of states. If the latter is a function of 
compact support, as in the case of harmonic crystals, then we obtain a 
representation of the partition function in terms of a power series of the 
type © P] . 

(b) Confluent hypergeometric function of the following type: 

*(M+l ; _x) = f;^( n + *) V (£>0), (9) 

n=0 n - \ / 

are expansions of the form (jSJ) since the terms [i n — ( n ^ e ) form a Haus- 
dorff sequence (see ©)■ 

(c) Hausdorff sequences can be constructed as follows ^H]- Consider a sequence 
{M"}o° °f ( rea l) numbers, and denote by A the forward difference operator: 
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Afi n = jti n +i — fi n - Then we have 
A k fi n = A x A x ■ ■ ■ x A fi n = ^ (-l) m QV„ +fc _ ro (fc = 0,l,2,...), (10) 

fc-timos m=0 

A is the identity operator, by definition. Now, suppose that there exists 
a positive constant M such that 

n 

(- + i) (p - 1} E 

It can be proved ^Hj that condition 1|11J1 is necessary and sufficient to 
represent the sequence {fin}^ as follows: fi n — J Q t n u(t) dt (see formula 
(@J), where u S £ p [0, 1]. Thus we can say that the set {fin}^ 1 , constrained 
by the condition ijTTl) . forms a Hausdorff sequence. 

The main purpose of the present paper consists in performing a Watson-type 
resummation of expansions of type (jSJ), where the set of numbers {Hn}^ is assumed 
to be a Hausdorff sequence generated by a function u(t) (see formulae (@J and dJ) 
which belongs to L^ 2+ ^[0, 1] (e > 0). In this case we can regard the sequence 
{/i n }g° as the restriction to the integers of a function fi{z) (z € C), which belongs 
to the Hardy space -ff 2 (C_i/2), and which is the unique Carlsonian interpolation 
[3] of the numbers {/i n }o°- We can thus perform the Watson-type resummation of 
expansion (JHJ) , and finally obtain another representation whose numerical handling 
is much more convenient and effective. 

The paper is organized as follows. In Section [21 we study the Carlsonian inter- 
polation of the Hausdorff moments {fi n }(f, and expand the function fi{iy— 1/2) in 
terms of the so-called Pollaczek functions. In Section [3] we perform a Watson-type 
resummation of expansion JHJ). In Section 01 we study an appropriate truncation 
procedure of the new representation obtained in Section |31 In Section we solve 
the Hausdorff moment problem by the use of the Pollaczek polynomials Ej , 
and show the connection between this problem and the Watson resummation of 
expansion JSJ. Finally, Sectionals devoted to numerical analysis and examples. 

2. Interpolation of Hausdorff moments and Hardy spaces 
We prove the following theorem. 

Theorem 3. Let the sequence {/j n }™ satisfy condition with p ^ 2 + e (e > 0). 
Then there exists a unique Carlsonian interpolation of the numbers fi n , denoted by 
fi(z) (z G C, fi(n) = fi n ), that satisfies the following conditions: 

(i) fi(z) is holomorphic in the half-plane Rez > —1/2, continuous at Rez = 

-1/2; ^ 

(ii) fi(z) belongs to L 2 {— oo,+oo) for any fixed value of Rez = x ^ —1/2; 

(iii) fi(z) tends uniformly to zero as z tends to infinity inside any fixed half-plane 
Rez ^ S > -1/2; 

Proof. If the sequence {/i n }o° satisfies condition with p ^ 2 + e (e > 0), then 

fln = [ t n u(t)dt, (12) 
Jo 



(-l)"- l A"-V, 



< M 



(n = 0,l,2,...;p>l). (11) 
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with u G L 2+e [0, 1]. Next, set t — e s in formula (|12(l . and obtain 

/>oo 

/!„=/ e~" s e- s u(e- s )ds (n = 0,l,2,...)- (13) 



n 



Therefore the numbers fi n can be regarded as the restriction to the integers of the 
following Laplace transform: 

/>OC 

H(z)= e- {z+1/2)s e- s/2 u(e- s )ds. (14) 



Indeed, one has fi(n) = [i n . By applying the Paley- Wiener theorem |10| to equality 
(|14f> . and recalling that the function e~ s / 2 u{e~ s ) belongs to L 2 [0,+oo), we can 
conclude that fi{z) belongs to the Hardy space H 2 (C_i/ 2 ), C_i/ 2 = {z G C, Rez > 
— 1/2} (sec Rcf. ^01)- We can thus state that fj,(z) is holomorphic in the half-plane 
Rez > —1/2, and tends uniformly to zero as z tends to infinity inside any fixed 
half-plane Rez S > —1/2. We can then apply the Carlson theorem 0, and say 
that (Jt(z) is the unique Carlsonian interpolation of the numbers ji n . Furthermore, 
in view of the fact that e~ s / 2 u(e~ s ) belongs to L 2 [0, +oo), then fi(— 1/2 + iy) 
belongs to L 2 (— oo, +oo), and, consequently, property (ii) holds true for any fixed 
value of Rez = x ^ —1/2. Finally, let us note that the function e~ s / 2 u(e~ s ) 
belongs to i 1 [0,+oo); in fact, \e~ s / 2 u(e~ s )\ ds = J Q \u(t)/y/t\ dt < oo since 
u G L 2+e [0, 1] (e > 0). Therefore, in view of the Riemann-Lebesgue theorem 
applied to representation i|14fl . it follows that the function n(— 1/2 + iy) (y G M.) is 
continuous, and thus property (i) is proved. □ 

Let us now introduce the following set of functions: 

*„( 2/ ) = -^rQ + ^^ l 1 /2) (2/) , ( i5) 

where T(-) denotes the Euler gamma function, and P^^ 2 \-) denote the Pollaczek 
polynomials P^\-), with A = 1/2 (see the appendix). These polynomials (in what 
follows the superscript A = 1/2 will be omitted) are orthonormal in (— oo, +oo) 
with weig ht function CHI 



w(y) = - 
ir 



1 * 2 



(16) 



Therefore the orthonormality condition reads: 

X) 

w(y)P n (y)P m (y) dy = <J n , m . (17) 



It can be proved that the functions {^niu)}^ form a complete basis in the space 
L 2 (— oo, +oo). Therefore the function //(— 1/2 + iy), which belongs to L 2 (— oo, +oo) 
(see Theorem|3J), can be expanded in terms of this basis. We can state the following 
proposition. 

Proposition 1. // the sequence {/Un}o° satisfies condition with p ^ 2 + e 
(e > 0), then 

1 \ °° 

2 + ^ ) = J2 ( 18 ) 



n=0 
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Figure 1: (A) Contour integration used for evaluating the coefficients c„ (see Theorem 
0J. (B) C + is the integration path of the integral in 12 H . 



iy Pn{y)dy. 



(19) 



which converges in the L 2 -norm. The coefficients c n are given by 

Proof. The sequence {Hn}"^ is a Hausdorff sequence satisfying condition 111|) with 
p 2 + e (e > 0), then /i(— 1/2 + iy) belongs to L 2 (— oo, +oo) (statement (ii) 
of Theorem El , and expansion l|18fl converges in the sense of the L 2 -norm; the 
coefficients c n are then obtained by the use of the orthonormality condition 1)170. □ 

The coefficients c n can be evaluated as follows. 

Theorem 4. The following equality holds true: 

.1)* 



= 2V^ 



OO 

£■ 

fc=0 



fc! 



(20) 



where P n {-) are the Pollaczek polynomials. 



Proof. Integral 1)19(1 can be evaluated by means of the method of contour integra- 
tion. Set in formula (|19|) : —1/2 + iy = z, and accordingly y = —i(z + 1/2). Then, 
performing an integration along the contour 7 shown in Fig. and taking into 
account the asymptotic behavior of the gamma function, we obtain 



7+ioo 



1 



1 



(i(z)T(—z)P n -i \ z + ^ j dz — J n(z)T(—z)P„ -i \^z + ^ J dz, 

(21) 

where C + is a path which encircles the real positive semi-axis of the z-plane in 
counterclockwise sense (see Fig. ^3). Then, using the theorem of residues we get 



c+ 



l x{z)T{-z)P n 



z + 



dz = 1\pn 



k=Q 



k\ 



f-kPn 



(22) 
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□ 

3. Watson resummation of a class of ^-transformed power series 

We prove the following theorem. 

Theorem 5. Expansion @ 7 where the terms {/^n}o° form a Hausdorff sequence 
satisfying condition (1 1 11 with p ^ 2 + e (e > 0), can be rewritten in the following 
form: 



fJ 

where 



(24) 



k=0 ' L V /_ 

P„(-) 6em<7 i/ie Pollaczek polynomials. The convergence of expansion H&gj l is uni- 
form on any compact subset of the real positive axis. 

Proof. We start by rewriting expansion (|SJ) in the following form: 

r l^nX n — r n n e na (a = \nx). (25) 

^ — ' rv. — ' n! 

n=0 n=0 

Next, using once again the theorem of residues, we rewrite the sum l|25(l as the 
following integral: 

E —T^e na = — / T(-z)^(z)e az dz, (26) 

where the path C + encircles the real positive semi-axis of the z-plane (see Fig. 
EJ3). Equality (|2*6*|l holds true since: 

(i) ju(z) is the Carlsonian interpolation of the moments {/i n }; 

(ii) for z = n (n = 0, 1, 2, . . .), one has /u(n) = //„, and the function T(—z) = 
r(— n.) is singular and has simple poles with residues (— l)"/n!. 

We can now close the contour C + as shown in Fig. We have, by exploiting the 
asymptotic behavior of the gamma function 

T(-z)^(z)e az dz = 0. (27) 

'7 

From l|27(l , and using the Stirling formula for the gamma function, 

/— 5+ioo 
r(-z)^i(z)e az dz 
-*-<<» (28) 

In the latter integral we use (I15|) and JTSJ. Then from (|26|l and (|28(l we have 
(-1)" 1 



E y -ih^ na = £ y r U - *v ) M * - 5 ) eQ( "" 1/2) * 

°° ' ' (29) 
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Using the formula 



iy T 



w 



cosh(7ry) 



from it follows that 

[ + °° ^ P n (y) 
/ 2^ Cn r 



cosh(7ry) 



(30) 



(31) 



where 5F denotes the Fourier integral operator. Interchanging integration and sum- 
mation, we have, from formulae pl - (|3Tl> : 



P n (y) 

cosh(7ry) 



o«xy dy = 



n=0 

-a/2 





1 


( _< s) 


cosh(a/2) 



20Fcosh(a/2) 



E 



Ctt, 2 



tanh 



(f) 



Substituting a — In x in (|32l) yields 

V2 



E ^-h»nx ri 



n=0 



n=0 



where 



k=0 



(-1)' 

fc! 



3C — 1 
X + 1 



-i\k + - 



(x > 0), 



(32) 



(33) 



(34) 



It remains to prove that the series at the r.h.s. of formula (|33(l converges uniformly 
on any compact subset of the positive real axis. Using the Schwarz inequality, 



E< 



x- 1 
x + 1 



Ei 



1/2 



E 



x + l 



1/2 



(35) 



The sum J^^Lo l M «| 2 can be shown to be finite by using the Parseval equality on 

2n 



expansion l(T%|l . The sum X)nLo 

' x-l \2 



x+l 



can be easily reduced to the series J2^=o V n > 

□ 



V = (f+r) ) which is uniformly convergent on any compact set y ^ yo < 1. 

4. Truncation of the resummed expansion 

We hereafter assume that only a finite number of Hausdorff moments fik (see 
formula Q34[ll are given, and, furthermore, we suppose that they can also be affected 

by noise, being typically round-off numerical errors. Accordingly, they will be 

(n) 

denoted by fi\ , f] denoting the order of magnitude of the numerical noise. Precisely, 



we state: 



1 1 fa) 



^ r\ (k = 0, 1, 2, . . . , k ; 77 > 0); (k + 1) is the number 

of Hausdorff moments which are supposed to be known. Next, we introduce the 
following finite sums: 



k=0 



k\ ^ k n 



i[k + 



1 



(36) 



s 
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With obvious notation we write: Un' 00 ^ = u n . Then, the following two auxiliary 
lemmas can be proved. 

Lemma 1. The following statements hold true: 



E| 

n=0 



,(0,oo) 



2 _ 1 
7T 



1 



dy = C (C = const); 



(ii) 



(iii) 



E| 

n=0 



-co; 



hm u<? ,fco) = M f °°) = u„ (n = 0, 1, 2, . . .); 

fco — >oo 



(iv) If mo(j],ko) is defined as 

{m 
m 6 N : l^n 7 '' 



n=0 



^ c 



(v) TTie sum 



lim 7710(77, fco) = +00; 

fco — >oo 



M(?' fe °) = ^|^' 7 ' feo) | (men), 

n=0 



satisfies the following properties: 

(a) It increases for increasing values of m; 

(b) the following relationship holds true: 



M (nM>) > 



1 



m — >-+oo 



(2m) 2fe ° (fc fixed). 



(fc !) 2 

Proof. The proof is given, with minor modifications, in Ref. g| 
Lemma 2. 77ie following equality holds true: 

mo ( J), fco) 2 

lim > 

bo — > oc ^ — ' 



I (, )fc0 ) _ u (0,oc) 



0. 



Proof. See Ref. g]. 

Next, we introduce the following notations 

00 



n\ x + 1 ' 

n=0 n=0 



x- 1 

X + 1 



y— m (r),fco) 
n=0 



X - 
fco 



fc=0 



fc! 



x- 1 

X + 1 



i[k + 



(37) 

(38) 
(39) 

(40) 
(41) 

(42) 



(43) 
□ 

(44) 
□ 

(45) 
(46) 
(47) 
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Then the following is true. 

Theorem 6. The following equality holds: 



lim 

ko — >oo 



g (v,ko) _ g(0,<X>) 



0. 



Proof. We have 



V2 



x + l 
V2 



mo(»),fco) 

E 



u^' ko) i n 



x - 1 



n oo 



i -y u ^ 



n=0 



x-l 

X+l 



X+l 



E 



,(0,oo) 



n=m a (r) ,k Q ) + l 
mo(ri,k ) 



E (^ hk0 



x- 1 
x + l 

x - I 
x + l 



Now, using the Schwarz inequality, 

'x-l 



E 



.(0,(30) 



n=m (r/,fe ) + l 



X+l 



1/2 



E 

i n=mo(rj,fco) + l 



,(0,oo) 



E 

<n=m (r], fco) + l 



x - 1 


2n\ 


x + l 





1/2 



Since 



and 



,(0,oo) 



< +oo, 



m (r],k ) > +oo, V 

/C(}— »OG ^ — * 

77— >0 "-=0 

oo 

^|(x-l)/(x + l)| 2 " <+oo forx>0 



n=0 



(see Lemma^l, the r.h.s. of formula l|5U|) tends to zero as ko — > +oo, 77 

m (r;,fc ) 



m o(»),feo) 

< E 

n=0 



x- 1 
x + l 

x-l 



' m (r/,fc ) 



From Lemma |21 it follows that 

m (ri,ko) 

lim y 



X+l 

1/2 



'm {r},ko) 

E 

n=0 



1/2 



X-l 



X+l 



u {r,,ko) _ u (0,oo) 



= 0, 
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while lim^^oo J2^o M 10 ~ l )/( x + l )\ 2n is finite for x > °- D 

77— v0 

5. Connection with the Hausdorff moment problem 

The classical Hausdorff moment problem can be formulated as follows ^ E| : 
Problem. Given a sequence of real numbers {fin}^, find a function u{t) such 
that 

fin = / t n u{t)dt (n = 0,1,2,...). (53) 
Jo 

This problem is ill-posed in the sense of Hadamard [Sj: Suppose, for instance, 
that we are looking for a solution in the space X = L 2 [0, 1], and assume that a 
solution in this space exists and is unique, but it does not depend continuously on 
the data. Further, in practical cases only a finite number of moments {fi n }o are 
known. We must then look for a solution in a finite-dimensional subspace X/v+i of 
X. Therefore, any function which is orthogonal to Xjv+i cannot be recovered: the 
solution is not unique. From the numerical point of view, we are led to the inversion 
of matrices which are severely ill-conditioned. We shall return on these questions 
later. For now we assume that a countable set of noiseless moments {/i n }o° are 
given, and prove the following theorem. 

Theorem 7. Suppose that the real sequence of Hausdorff moments satisfy 

condition ll 11\ with p 2 + e (e > 0). Then the function u(t) can be represented by 
the following expansion, which converges in the L 2 -norm: 



00 (—l) k 

U{t) = u n$n{t), U n = \/2^ fc , flkPn 

n=0 fc=0 



(54) 



where 

*„(*) = i n V2e- 1 L n {2t), (55) 
and -P n (-) and L n (-) are the Pollaczek and the Laguerre polynomials, respectively. 

Proof. In formula (|14f> set z = —1/2+iy, Recalling that the support of the function 

f(s) = e - s / 2 u{e- s ) belongs to R+, 

^(~\ + ly )^J + e ~ tVSe ~ s/2u ^ S ) ds = 3 r {e" s/2 u(e" s )} , (56) 

where 3 denotes the Fourier transform operator. Let us now return to the expansion 
(I18[) and to formula (|15() . which gives the expression of the functions ^f n (y). In 
particular, in the integral representation of the Euler gamma function T(l/2 + iy): 
i.e, T(l/2 + iy) = f+°° e^^ 1 / 2 ) dt, we set t = e~ s : 

v(^+iy\=J e- e ~ s e- s/2 e- lsy ds = ?{e- e ~ s e- s / 2 y (57) 
Since the function e~ e e~ s l 2 belongs to the Schwartz space S^, 



ds 



^-e^ s e- s ' 2 



(58) 
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therefore, from equality (|56|l . and recalling once again expansion (|18J) . we obtain 



Reverting to the variable t = e~ s , we have: 

oo 

«(*) = E 



- S /2 



V n Vdt 



Vte~ 



(59) 



(60) 



71=0 



where u„ = c n /^/2n, and the functions 3> n (£) are given by 

$ n (t) = V^Pn f*t^ [Vte"*] = V2 l "e- t L„(2t), (61) 

where the £„(•) are the Laguerre polynomials. Note that the functions $> n (t) form 
a complete basis in L 2 (0, oo), and that the convergence of expansion I^U|) is in the 
sense of the L 2 -norm JS]. □ 

We now analyze the truncation of expansion (|6U[) . We introduce the approxima- 
tion 

™o(?),fco) 



E 



$n(*), 



(62) 



where the coefficients u 
formula (@DJ. 

Theorem 8. 



Proof. 



are given by formula l|36|) . and mo (77, &o) is defined by 



lim 

fco — >oo 



L 2 [0,+oo) 



u (»J.fco) _ u 



L 2 [0, + oo) 



E 



E 



,(»7.feo) 



n=m (r),feo) + l n=0 

The statement of the theorem follows from Lemmas ^ and |21 



(63) 

(64) 
□ 



As already remarked, the Hausdorff moment problem, formulated as above, is 
severely ill-posed. In principle one could use regularization procedures [SJ I12| . 
among which the Tikhonov's or Tikhonov-based methods are the most popular 
\T%\. Each of these procedures consists in restricting the class of admissible 
solutions to a compact subspace of the solution space (for instance, a subspace 
of X = L 2 [0, 1]), by introducing suitable bounds on the solutions. However, some 
problems remain, and in particular the determination of the so-called regularization 
parameter, whose optimal choice requires a precise knowledge of the majorizations 
on the solutions and on the noise affecting the data. 

The method presented above does not make use of any a-priori knowledge on 
the solution and on the data. In several cases the truncation given by formula l|62|) 
can be easily determined by the properties of the truncation number mo (77, fco), 
illustrated by Lemma ^ and, in particular, by the statements (iv) and (v) (for 
numerical examples, see the next section). Then, the statement of Theorem [S] 
guarantees the convergence of approximation l|62|) to the solution, as the number 
of data increases to infinity, and the noise tends to zero (see H63[l). If the number 
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of data is too small, or the noise is too large, or, finally, if the function u to be 
determined is irregular (i.e., presents discontinuities of various types), the sum 
Mm' k °' > (see Q431 in LemmaQJ can present no plateau (see next section), and the 
method cannot be used. However, this negative result still provides information: the 
continuity which could at best be restored with classical regularization procedures 
remains extremely weak. 

Returning to the Hausdorff-transformed power series of type JBJl, it is worth 
remarking on the following fact: the method of resummation which we have pre- 
sented is affected by the same type of ill-posedness illustrated above in connection 
with the solution of the Hausdorff moment problem. More precisely, we face the 
ill-posedness connected to the reconstruction of the function /i(iy — 1/2) from the 
sequence of Hausdorff moments {fik}o ■ We can thus advance the following critical 
remark: the method of resummation presented transforms the type of pathology 
affecting expansion JHJ), i.e., slow convergence, into another type of pathology, i.e., 
ill-posedness. However, the ill-posedness of the problem, cured by the truncation 
procedure presented above, has, at least in the case of regular Hausdorff transfor- 
mation, much milder effects than the slow convergence pathology on the actual goal 
to be achieved: the numerical evaluation of functions of type (JSJ). In other words, 
the regularization of the ill-posedness cures the drawbacks of the slow-convergence. 

6. Numerical analysis 



Return to formulae lf2*3|) . J2U, an d set x = 1/r, (r e [l,+oo)). Then 




Next, we take for c the value r G [1, +oo); then the rightmost expansion in (|65fl and 
formula (|66|l reproduce once again the original series X^^Lo n\ ^" (l/ r )"j since 
Pq(') = 1- The same type of result can be obtained when the number of moments 
fib is finite. We can thus conclude that the expansion that we have proposed can 
always be reduced to the standard one in the interval (0, 1], where the original series 
converges rapidly. We shall now show that the resummed expansion converges much 
more rapidly than the original one for x ^> 1. First we discuss in detail how to 
manage this new expansion numerically, and, in particular, how the truncation 
number rao(»),fco) can be determined. 

The problem of evaluating mo(?7, fco) is intimately related to the price that must 
be paid for coping with the ill-posedness of the analytic continuation involved in the 
reconstruction of the function fi(iy — 1 /2) from the sequence of Hausdorff moments 
{/ifcj-Q . As remarked in Section[Sl most regularization procedures generally require 
a priori bounds on the solution and on the data. Then one is led to introduce 
in the regularizing algorithm a truncation, or a filtering, which depends on the a 
priori information on the solution that one is supposed to have. Our procedure 
does not require any prior knowledge, and the truncation number 777.0(77, ko) (see 
Section 0J) can be determined by analyzing the sum versus m (see |@2J). 
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From statement (v) of Lemma and from formula (|41|l it follows that, if fco is 
sufficiently large and r\ is sufficiently small, then Mm ,k °^ presents a plateau, and, 
after that, it starts growing as a power of (2fc ) (see and Figs. |23, [3p and 
^\.) . The extension of the plateau varies, and increases as the number of moments 
lik becomes larger (see Fig. EJ3). Analogous results are obtained even in the case 
?/ 7^ (see Fig. 0K). A simple algorithm for the automatic derivation of mo can be 
implemented along the following lines. We start by observing that the knowledge 
of the asymptotic behavior of Mm' k °\ for large m, allows us to restrict the range 
of mo by defining an upper limit m a (mo < m a ), which represents the value of m 
where approximately the asymptotic behavior sets in; in practice, m a is set as the 
value of m where Mm' k °^ starts being close enough to its asymptotic behavior. The 
candidate plateaux are then located by selecting the extended intervals of m < m a 
where the modulus of the first numerical derivative of Mm' k °^ is sufficiently small. 
Finally, mo is chosen as the largest value of m belonging to the interval which is 
closest, but inferior, to m Q . ft should be noticed that the choice of mo within 
the plateau is not critical for the accuracy of the final result. For the sake of 
completeness, it should also be mentioned that the erratic behavior of the noise 
can produce very short plateaux located between the true value of mo and before 
Mm' fe °' starts following its asymptotic behavior (i.e., for m = m Q ). In this case 
our procedure could fail to recover the correct value of mo; this drawback has been 
solved heuristically by simply rejecting plateaux shorter than a given threshold 
length. 

We test our method by comparing different evaluations of the confluent hyper- 
geometric function $(1,£+ 1;— x) (see formula ©)■ The latter can be directly 
evaluated as the Laplace transform of the following function: 



On the other hand, one can evaluate $(1,£ + l;—x) by means of the standard 
expansion at the r.h.s. of formula jnj (truncated at a certain n — no), and finally 
these results can be compared with those obtained by the Watson resummation 
method (in particular, see H4tj[) ). The results are illustrated in Figs. 0E1 

In Fig. [2K the plots of $(1, 2; —x), computed by using the standard formula © 
with different values of no, are shown. It is evident how the deviation from the true 
function rapidly explodes as x increases. Moreover, even using more moments fi n , 
i.e., increasing no, the situation does not get better significantly. Figure I2J3 shows 
the sum Mm ,k °^ versus m for various values of fco (&o = 5, 10, 15,20). It can be 
seen the presence of the plateaux, whose length increases as the number of moments 
Hk used in the computation, i.e., fco, increases (see ()41|l in Lemma [Q. This figure 
shows how the truncation number mo(0, fco) can be determined. The comparison 
among the true <I>(1, 2; —x), computed analytically as the Laplace transform of the 
function u(t) in I|tj7fl (solid line), the evaluation obtained by the truncated standard 
expansion © (dashed line), and truncated resummed expansion i|46|) (filled dots), 
is shown. In this case, fco = no = 20. From the inspection of M™' 20 '' in Fig. 
03, it is recovered that the plateau approximately ranges from m = 8 through 
m = 38. The computation shown in Fig. 03 was made with mo(0,20) = 30. It 
is evident how the accuracy of the computation increases considerably when the 




te [o,i], 
t e (i,+oo). 



(67) 
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Figure 2: Computation of the confluent hypergeometric function + 1; — a;) with 

1 = 1. (A) Computations of $(1, 2; —x) by means of expansion truncated at n — no, 
for different values of no. The number near each plot indicates the value of no used for 
the computation. The solid line represents the actual function $(1,2;— x). (B) Jfi '' ' 
versus m for fco = 5, 10, 15, 20 (see 14211 ). (fco + 1) is the number of moments a k used for 
the computation. (C) $(1, 2; —x) computed by means of the resummed expansion 1461 
(dots) by using 21 moments u k (i.e., fco = 20). mo(0, fco) = 30 has been set according to 
the analysis of the function Afi 0,20 ' in (B). The solid line represents the actual function 



$(1,2; 



The dashed line shows the computation made by using expansion @ trun- 



cated at no = 20. (D) Plot of the relative error of reconstruction |($ a 



$truc)/$t 



in the range x (E [1 : 100]; no = fco = 20. The solid line shows the relative error computed 
by using expansion @. The dashed lines indicate the relative error of the computations 
made with the resummed expansion 11461 . with mo(0, 20) = 6, 30, 45. 



resummed expansion is used. This fact is made even more clear in Fig. [2)3, where 
the relative error of computation |A$/$| = |($ ap pr — 3>truc)/$truc| over the range 
x G [1 : 100] is shown. The error made by using the standard expansion (solid 
line) diverges, whereas the error made by using the resummed expansion (dashed 
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line) remains quite limited. This latter is displayed even for the cases too — 6 and 
mo = 45, which represent two values just outside the plateau. In these cases the 
error increases, though in different ways, with respect to the correct value mo = 30. 
For any other value of mo within the plateau, that is to £ [8,38], the error of 
computation does not change significantly from that at too = 30. This fact is 
evidentiated in Fig. EK, where the root mean square error of the computation of 
$(1,2;— x), integrated over the interval x S [0,20], is shown for ko = 20, and Too 
varying in the range [0 : 50]. As long as mo remains into the plateau the error is 
small, whereas when too increases far beyond the upper limit of the plateau, the 
error rapidly grows. In order to show the interplay between ko and too, in Fig. 
03 it is shown the relative error for different values of ko {ko = 0, 4, 8, 10, 14), and 
too = 10 kept fixed. Even in this case the error is high when too does not lie 
within the plateau of Afi°' fco) (see also Fig. Hp). In Fig. Ep the sums M^ 2a) have 
been plotted for various values of £ (£ = 1, 2, 3, 4, 5, 8, 15, 20), while in Fig. the 
corresponding computations of $(1,£+ l;—x) (£ = 5,8,15,20), made by means 
of the resummed expansion 1)46(1 with ko = 20 and Too (0,20) = 10 are compared 
with the true functions. From the analysis of Ml"' 20 ' in Fig. [3p it can be seen 
that the value to o (0, 20) = 10 lies outside the plateau resulting for £ = 15 and 
£ = 20; correspondingly, in Fig. the computation of $(1, £ + 1; —x) with £ = 15 
and £ = 20 clearly deviate from the corresponding true functions. Finally, Fig. 0] 
illustrates the analysis in the case of noisy moments //j^ • To obtain the , the 
moments fik have been noised by adding white noise uniformly distributed in the 
interval [—77,77]. Figure 0J\ shows the sum Mt W) for various values of the noise 
parameter 77 (77 = 10" 1 , 10~ 2 , 10~ 3 , 10~ 4 , 10~ 5 , 10~ 6 ; £ = 1). These plots show how 
the plateau gets shorter as the noise level increases. In Fig. some examples 
of computation of $(1,2;— x) by using the standard expansion for various values 
of r\ are shown, while in Fig. ^JP the computations made by using the resummed 
expansion are given. This panel shows that a significant deviation from the actual 
function arises only in the case of quite noisy moments (77 ~ lO^ 1 ), and remains 
quite acceptable for small levels of noise. The difference of accuracy achieved by 
using the two types of expansion are made evident in Fig. where the root mean 
square error, in the range x € [0, 10], is shown as a function of r\. The dots indicate 
the error referred to the standard expansion, while the crosses indicate that referred 
to the resummed expansion. Over a wide range of noise level, the error made by 
using the resummed expansion remains orders of magnitude smaller than that made 
by using the standard expansion. 



The main properties of the Pollaczek polynomials P,[ (y) are briefly summa- 
rized: 

(1) In terms of the hypergeometric series 0, 



Appendix 
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(A.2a) 
(A.2b) 
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Figure 3: (A) Root mean square error of the computation made with the resummed 
expansion 1461 with respect to the actual function $(1,2; —x) in the range x £ [0 : 20], 
versus mo(0, fco); fco = 20. (B) Relative error of computation of $(1, 2; —x) by using the 
resummed expansion 114611 in the range x £ [1 : 1000]. mo = 10 is kept fixed, while fco 
varies: fco = 0, 4, 8, 10, 14. (C) Computation of the confluent hypergeometric function 
$(l,f + 1; — x) for different values of I. In this panel the plots of Mi 0,20 ' (i.e., fco = 20) 
versus m are shown. Each plot has been computed by using the moments fik = 
with £ = 1,2, 3, 4, 5, 8, 15, 20. (D) Plots of $(1, l+l; —x) computed by using the resummed 
expansion 1461 (dots), for various values of £: £ = 5, 8, 15, 20; fco = 20 and mo(0, 20) = 10 
have been kept fixed. The lines represent the actual $(!,£+ 1; —as). 



(3) The generating function is given by: 

oo 

z n Pi 1 ' 2 \y) = (1 - izfy~^ 2 \l + izf-^-W (|*| < 1). (A.3) 

n=0 
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Figure 4: Computation of $(1, 2; — x) by using noisy moments (i£ . The moments \ik 
have been noised by adding white noise uniformly distributed in the interval [—77,77]. (A) 
Mm' 20 ^ versus m computed for n ranging from n = 10 _1 through r\ = 10~ 6 with step 
10 _1 . The rightmost solid line indicates the noiseless Mi"' 20 '. (B) Comparison between 
the actual function $(1, 2; —a;) (solid line) and the computations made by using expansion 
©, with 77 = 10" 2 , 1CT 3 , 1CT 4 . The dashed line labelled by "0" has been computed by 
using the noiseless moments fik- The number of moments used is (no + 1) = 21. (C) 
Comparison between the actual function $(1, 2; —x) (solid line) and the computations 
made by using the resummed expansion 1461 : fco = 20. The filled dots indicate the 
computation made with 77 = 10~ 3 and mo(10~ 3 , 20) = 5 (see also panel (A)). The dashed 
lines represent the computations made with 77 — 10~ 2 , mo(10 _2 ,20) = 4, and r\ — 10 -1 , 
mo(10~ 1 ,20) = 2. (D) Root mean square error of the computation of $(1,2;— x) versus 
r] in the range x G [0 : 10]; fco = no = 20. The dots indicate the error made by using 
expansion ©, while the crosses represent the error made by using expansion 14611 . The 
truncation numbers are: m (10" 6 ,20) = 30, m (10" 5 ,20) = 15, m (10~ 4 ,20) = 10, 
m (10" 3 , 20) = 5, m (10~ 2 , 20) = 4, mo(10 _1 ,20) = 2. 
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